Subdiffusion and cage effect in a sheared granular material 
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We investigate experimentally the diffusion properties of a bidimensional bidisperse dry granular 
material under quasistatic cyclic shear. The comparison of these properties with results obtained 
both in computer simulations of hard spheres systems and Lenard- Jones liquids and experiments 
on colloidal systems near the glass transition demonstrates a strong analogy between the behaviour 
of granular matter and these systems. More specifically, we study in detail the cage dynamics 
responsible for the subdiffusion in the slow relaxation regime, and obtain the values of relevant time 
and length scales. 
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Glass forming systems exhibit many intriguing prop- 
erties, and their study has generated a large amount of 
theoretical as well as experimental work. One of the main 
features of their dynamics is what is usually called the 
cage effect, which accounts for the different relaxation 
mechanisms^, 0: at short times, any given particle is 
trapped in a confined area by its neighbours, which form 
the so called effective cage, leading to a slow dynamics; at 
sufficiently long times, the particle has managed to leave 
its cage, so that it is able to diffuse through the sample by 
successive cage changes, resulting in a faster relaxation. 
These mechanisms define the (3 and a regimes. 

Many experiments and simulations have been per- 
formed to study this scenario. In particular, the under- 
standing of the nature of cages require microscopic in- 
formations which have been essentially obtained in com- 
puter simulations of hard spheres systems and Lennard- 
Jones liquids (e.g. 0, S IH IS ) ■ A suitable way to ex- 
tract these informations in laboratory experiments con- 
sists in using systems undergoing a glass transition com- 
posed of sufficiently large particles so that it is possible to 
follow them through direct observation. The main exam- 
ple of this method is the breakthrough study of colloidal 
particles near the glass transition by confocal microscopy, 
realized by Weeks et at 0, 0, El , who first observed ex- 
perimentally the cage effect in real space. 

Beside, especially since the crucial experiments of the 
Chicago group 0, 0] , it is widely supposed that dense 
granular matter could be considered as an analog of 
glassy systems, albeit a rather special one, since it is 
athermal |l3| . Granular systems also undergo a jamming 
transition which shares many properties with the glass 
transition, arising the possibility of a unified description 
\lA UM ■ Then, granular materials could represent a 
simple way to perform accurate measurements and un- 
derstand the nature of cages (since grains can be rela- 
tively lar ge a nd then quite easy to follow through direct 
imaging, [rjj) provided one checks that the analogy no- 
ticed at the scale of the sample is confirmed by a precise 
study of the diffusion properties at the grain scale. 



In this paper, we show that a very simple system 
such as a bidimensional bidisperse dry granular material 
submitted to a quasistatic cyclic shear indeed behaves 
the same way as glassy systems in the sense that its 
diffusion properties evaluated with the same tools as in 
numerical studies of glasses 0, 0- B behave identically 
(except for the initial thermal regime which obviously 
does not exist in a granular material). As a result, 
we also characterize the cage effect and the dynamical 
heterogeneities. 




FIG. 1: Experimental setup (left: picture; right: scheme). 

The experimental setup is as follows (figHJ : a bidimen- 
sional, bidisperse granular material, composed of about 
6.000 metallic cylinders of diameter 4 and 5mm in equal 
proportions, is sheared quasistatically (7 = jd/^P/p = 
10 -5 ; see |l8|) in an horizontal deformable parallelogram. 
The shear is periodic, with a shear rate 7 = 1.5°.s _1 and 
an amplitude 6 max = ±10°. The volume accessible to 
the grains is maintained constant by imposing the height 
of the parallelogram, so that the volume fraction is a 
constant (<f> = 0.86). We follow a sample of 500 of the 
grains with a CCD camera which takes a picture of the 
material each time the system is back to its initial posi- 
tion (0 = 0). The unit of time is then one cycle, a whole 
experiment lasting 10.000 cycles. The unit of length is 
chosen to be the mean particle diameter d. 

Fig El 3 shows a typical trajectory: the particle spends 
most of its time confined in a well defined area and some- 
times escapes during rare and brief events. In the follow- 
ing, we will refer to this behavior as " cage effect" . 
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FIG. 2: (a)some tracers trajectories (b)a typical trajectory; 
the circle indicates the particle size (c)grey: the same tra- 
jectory; black:part of the trajectory, showing the existence of 
cages at a scale of order 0.3 diameter. 



In order to precise the nature of this motion, we first 
study the statistical properties of the displacements dur- 
ing r time steps: AX{t) = X(t + r) - X(t). The pdf of 
these displacements are presented on figEt-, for r =1, 10, 
100, 1000. They are caracteristic of interminent dynam- 
ics, with fat tails compared with the gaussian case, and 
then can be interpreted as the signature of the cage ef- 
fect, which sidesteps the relative proportion of small and 
large events. 
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FIG. 3: (a) pdf of AX(r)/a for r 

1(.), 10(*), 100(o), 1000(+); the solid line is the gaussian 
distribution (b) ct(t) = (Ar 2 (r)); dotted lines show the 
slopes 1/4 and 1/2; dashed lines indicate the position of the 
crossover which determines r* and t* (inset: non-gaussian 
parameter a(r)). 



The root mean square displacement ct(t) = \/ Ar 2 (r) 
presents two regimes(fig|3|D) : at short times, the dynam- 
ics is subdiffusive (logarithmic slope 1/4), which can be 



interpreted as the result of the trapping of the particles 
in cages during the /3-relaxation, while it becomes diffu- 
sive (logarithmic slope 1/2) at long times, when particles 
have succeeded in escaping from the cages, leading to the 
a-relaxation. The crossover between the two regimes can 
then be considered as a measure of the cage size r* and 
cage lifetime t* (see fig Et*) ■ Here, we find r* ~ 0.3 and 
t* ~ 300. Note that this value of t* lies at the end of 
the range where the non-gaussian parameter a, defined 
as a(r) = (AX 4 )/(3(A1 2 ) 2 ) - 1, is maximum (inset of 
fig ED 3 )- This can be understood this way: on timescales 
much shorter than t*, as time grows, the particles explore 
their cage, being statistically more and more influenced 
by cage borders; when time approaches t*, they have ex- 
plored the whole cage and have constructed their statis- 
tics, so that a does not change significantly; on timescales 
longer than t* , they diffuse from cage to cage, and are 
then less and less influenced by the effect of trapping, so 
that a decreases. 

Note that the value of r* is smaller than the one which 
could be infered by direct reading of particles trajectories 
(figHb)- To make this clearer, we then choose to plot only 
a part of the trajectory. The result is shown on fig 0:. It 
appears that a region that we would visually describe as 
one cage is in fact the superposition of subcages, each of 
them having a size of about 0.3, i.e. r*. 

Following Doliwa & Heuer 0, 0] , we now turn our at- 
tention to the conditional probability P(xi2|rrji; r) (resp 
P(j/i2|roi; t)), which represents the probability distribu- 
tion of the projection X12 (resp. yu) of the motion during 
a time interval r along (resp. orthogonaly to) the direc- 
tion of the motion during the previous time interval, un- 
der the condition that the length of the motion during the 
previous interval has the value roi . Results are shown in 
figEl This quantity contains many informations 0: (i) at 
a given roi , the distributions are symmetric around their 
mean value; (ii) the distributions are not gaussian; (iii) 
the mean value of j/12 is 0, while the mean value of X12 is 
always negative; (iv) if we focus on the evolution of (X12} 
with roi, we observe two regions: for roi < 0.3, there is 
a linear relation between (X12) and roi ((^12) = c(r)roi, 
c(t) < 0), wherehas for roi > 0.3, (£12} is a constant. 
Note that the value of r* measured using property (iv), 
i.e. by the localisation of the crossover between the linear 
and constant evolutions of (0:12) with roi, is the same as 
the one extracted from the measurement of c(r). 

The slope \c(t)\ decreases with r(figEt,d), and is ap- 
proximately related to the logarithmic slope 8(t) of <t(t) 
by Hi: 

S est (T) =0.5 + ln[l + c(t)]/ In 4 

FigEJi shows both S(t) computed directly from <r(r) and 
the one computed using this formula. One sees that the 
behavior is well reproduced, despite a little offset which 
might be due to the approximations made in the calcu- 
lation of 5 es t- 
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FIG. 4: (a) and (b) Conditionnal probabilities (in colorscale) 
P(xi2|roi; t) and P(j/i2|roi; t) (see text); the white traces are 
the mean values (£12) and (1/12} (c) (2:12) for different values of 
t ( from bottom to top: r = 100, 300, 500) (d) plot of S(t)(.) 
and (5 es t(r)(A) (inset: slope c(r)). 

To go further in the interpretation of these distribu- 
tions, we extract their widths an and o~±_ . Their evolu- 
tion with roi is shown on fig Eli, for two different values 
of the time interval r. First, we note the increase of a / / 
with 7*01. This means that large steps are more likely 
for particles which moved farther during the previous in- 
terval. It is not the case in a purely diffusive process, 
since large events are statistical effects, with an occurancc 
which is not related to the length of the previous step, 
making all the particles to be equivalent. Here, since par- 
ticles which move farther are the ones which were already 
making large steps, this shows the existence of a popu- 
lation of fast particles, which is a typical feature of glass 
forming systems, as pointed out for example in 0, 0, - 
Second, we see that for short time intervals r, the in- 
crease of a J 1 is larger than the one of <r±. This reflects 
the anisotropy of the motion, like stringlike cooperation 
observed numericaly by Donati et al. [3. Both effects 
concern movements on short timescales, since they tend 
to disappear as we increase the time interval r (fig'EJi). 

At this point, we can give a partial conclusion about 
the diffusion properties of the system. The fact that 
the linear regime of (xi2)(roi) ends at the same value 
roi ~ 0.3 for any values of r < t* indicates that, in this 
regime, it is necessary to describe the system as driven 
by spatial constraints. For displacements smaller than r*, 
the larger a step the more anticorrelated is the following 
step, as shown by the negative values of c(t), which re- 
flects the systematic backdragging effect experienced by 
the particle in its cage. For displacements longer than 




FIG. 5: (a) Widths of the distribution of X12 (°7/) an d 1/12 
(er±) versus r i for r = 10 and r = 500 (b) [F 2 (t)] 2 , F 2 {2t) 
and Fs(t,t) (see text); at short times, F3(t,t) = [^(t)] 2 , and 
at long times, F 3 (t,t) = F 2 (2t). 

r* , a cage rearrangement has occured, so that the anti- 
correlation does not increase any more. The constancy 
of (£12) at this saturation value shows a memory of the 
fact that part of the trajectory was made in a cage.|c(r)| 
decreases with r, i.e. cage effect becomes weaker, which 
shows that cages relax and adapt to the new position of 
the enclosed particles. The fact that the typical relax- 
ation time is t* is justified by the need for the particles 
forming the cage to escape their own ones in order to 
adapt |3 • On timescales longer than t* , these effects dis- 
appear, which is the result of the increasing number of 
particles which have undergone a rearrangement. 

We now discuss the dynamical heterogeneities by con- 
sidering multitime correlation functions. Let us de- 
fine the three quantities F2(t) = (cos (q.f 01)) , F%(2t) = 
(cos(q.f 02 )) and F 3 (t,t) = (cos (q.f 12) cos (q.f 01)} , where 
the vector is the displacement observed between the 
times i x t and j x t: fij = f(j x t) — f(i x t). It 
has been shown that one can decide wether the dynamic 
is heterogeneous or homogeneous by comparing F^(t, t) 
with F 2 (2t) and [F 2 (t)] 2 respectively EllliJ. This can 
be understood by considering the definitions of homoge- 
neous and heterogeneous dynamics jl^ : in the purely 
homogeneous case, the movements during two subse- 
quent time intervals along a given direction q are not 
correlated in length, whereas in the purely heteroge- 
neous case, they are only correlated in length. Then, 
in the homogeneous case, since a cosine is sign inde- 
pendant, both terms of the product are uncorrelated, so 
that one can factorize Fs(t,t) and obtain [^(f)] 2 . In the 
heterogeneous case, one can replace cos(q.fi2)cos(q.foi) 
by cos(q.fi2 + q.roi) + sin(q.fx2)sin(q.roi). As in this 
case the signs are not correlated, the mean of the sec- 
ond term must be 0, so that we are left with F 2 (2t). 
Note that for a random walk, both equalities are ful- 
filled, so that F 2 (2t) — [F 2 (t)} 2 , i.e. the relaxation is 
exponential. These functions are presented in fig [5Jd, 
where we have chosen q = 27r/r*. One sees that at 
short times (t < 10) the dynamics is mainly homogeneous 
(F 3 (t,t) = [F 2 (t)] 2 ), and then slowly evolves toward an 
heterogeneous dynamics as time grows. 
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To better characterize the cooperation in the system, 
we use a convenient tool proposed by Hurley et al. 0, 
based on relaxation times. For a particle i, the relaxation 
time Ti(r) is defined as the time needed by the particle 
to reach a given distance r for the first time. The distri- 
bution of these relaxation times is shown in the inset of 
fig El 1 : f° r r = 0.3. The idea of Hurley et al. is the fol- 
lowing: if one defines Tij(r) as the mean relaxation time 
of the particles contained in a circle of radius / centered 
on particle i, then the study of the difference between 
Ti.i(r) and T av (where T av is the mean relaxation time 
calculated over all the particles) should give some infor- 
mations about the typical length L over which coopera- 
tive effects take place. The simplest quantity to compute 
is then the second moment (in the notations of 0): 



m 2 (l) 



{(T it i - T av f) 

<(t m - T av y) 




FIG. 6: (a) second moment rri2(l) of the relaxation time 
distribution, for different values of the cutoff distance r (0.1(.), 
0.3(*), 0.5(+)); the dependence of these curves on r is not 
monotonic (inset: relaxation time distribution for r = 0.3) 
(b) characteristic length L; it has a maximum L* for r ~ r*. 

7712(0 is plotted on figHO for different values of r. One 
clearly sees that the typical length on which 7712 decreases 
has a maximum around r ~ 0.3. To quantify this, we plot 
L (defined as the integral over I oi m 2 ) versus r and ob- 
tain the curve of fig^- Two important informations can 
be deduced: first, as we already noticed on the curves of 
m%(f), L as a maximum for r ~ 0.3, i.e. the typical cage 
size. This means that cage rearrangements are phenom- 
ena which imply more cooperativity than the dynamics 
at other scales (a complete discussion is given in ref @). 
Second, we obtain a value for this typical length L which 
is, at its maximum, L* ~ 7.5 particle diameters. We then 
see that cage rearrangements are highly cooperative phe- 
nomena. This, added to the small value of r* , shows that 
the picture of a particle escaping from its nearest neigh- 
bours is definitely not adapted. Instead, these events are 
subtle and complicated rearrangements, involving a large 
number of particles. One remaining open question is the 
nature of the slow phenomenon to which the cage dy- 
namics participate, since, contrary to the experiment of 
Pouliquen et al. [ljj on granular compaction, the present 



one is not submitted to gravity and has a constant vol- 
ume. Such aquestion as already been raised by Kabla & 
Debregeas |21| in the case of a gently vibrated pile. 

The first conclusion of this experimental work is that 
scenari imagined to describe glass forming systems ap- 
ply to the peculiar case of dense granular materials, giv- 
ing a precise sense to the strong analogy between these 
two fields and contorting the application of theoretical 
ideas inherited from statisical physics of glassy systems 
to granular matter. The second important result is that 
the peculiar dynamics observed in glassy materials still 
exists in a system to which an agitation far from a ther- 
mal noise is provided. Also, this experiment shows that 
a simple granular system as the one described in this pa- 
per may be an efficient laboratory model to look for the 
microscopic phenomena involved in cage rearrangements, 
which are still poorly understood, but of capital interest 
in the understanding of jamming and glass transitions. 
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